A model of color confinement* 



Daniel Zwanziger^ 
^Physics Department, New York University, New York, NY 10003, USA 

A simple model is presented that describes the free energy W{J) of QCD coupled to an external 
current that is a single plane wave, J{x) — H cos{k ■ x). The model satisfies a bound obtained 
previously on W{J) that comes from the Gribov horizon. If one uses this model to fit recent lattice 
data — which give for the gluon propagator D{k) a non-zero value, -D(O) 7^ 0, at fe = — the data 
favor a non-analyticity in W{J). 
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I. INTRODUCTION 

Recent numerical studies on large lattices of the gluon 
propagator D{k) in Landau gauge in 3 and 4 Euclidean 
dimensions, reviewed recently in [l| , yield finite values for 
-D(O) [3 - 0, in apparent disagreement with the the- 
oretical expectation that D{0) = 0, originally obtained 
by Gribov ^|, and argued in Q. Upon reviewing the 
argument Q which leads to D{0) = 0, one hypothesis 
stands out which should perhaps be dropped in view of 
the apparent disagreement. This is the hypothesis that 
the free energy W{J) in the presence of sources J is an- 
alytic in J. This is an important point because a non- 
analyticity in the free energy is characteristic of a change 
of phase. 

The free energy W{J) enters the picture because it is 
the generating function of the connected gluon correla- 
tors. In particular the gluon propagator D^.y, is a second 
derivative of W{J) at J = 0, 



D 



d^W{J) 



x,y 



dJxdJy 



(1) 



Here a condensed index notation is used, where Jx repre- 
sents J^{x), /i is a Lorentz index, and 6 is a color index, 
and we write 

(J, A)^Y1 = / J>Ki^)- (2) 

The free energy W{J) in the presence of sources J is 
given by 

expW^(J) = (exp(J,A)) 

= J dA p(A)exp(J, A), (3) 



where p{A) is a positive, normalized probability distri- 
bution. The integral over A is effected in Landau gauge 
c^M^M ~ ^- "^^^ domain of integration is restricted to the 
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Gribov region D,, a region in ^-space where the Faddeev- 
Popov operator is non-negative, M{A) = —d^D^{A) > 
0. 

The model will be defined for the special case where 
the source is given by 



J^(x-)=if^os(fc.x), 



(4) 



so the free energy 



expWkiH) = ( exp[^i7;:cos(fc.a:)4(a;)] ), (5) 



depends only on the parameters fc^ and Hj[. Because 
Afj^(x) is transverse, only the transverse part oif H is 
operative, and we impose 



0. 



(6) 



By analogy with spin models, may be interpreted as 
the strength of an external magnetic field, with a color 
index b, which is modulated by a plane wave cos(fc • x). 

A rigorous bound for Wk (H) on a finite lattice is given 
in Q. One can easily show that in the limit of large 
lattice volume V, and in the continuum limit, this implies 
the Lorentz-invariant continuum bound in D Euclidean 
dimensions. 



(7) 



where 



invariant norm of H^,. 



''^'^ is the color- and Lorentz- 
In our notation the vector po- 



tential is given by A{x) = g A^'^'^*^ (x) , and has engineering 
dimension in mass units [^(a;)] = 1 in all Euclidean di- 
mension D, while [H] = D — 1.] 

This bound yields in the zero-momentum limit 



wo{H) ^ lim WkiH) 



0. 



(8) 



As discussed in [oj, this states that the system does not 
respond to a constant external color-magnetic field no 
matter how strong. In this precise sense, the color degree 
of freedom is absent, and color is confined. 
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II. THE MODEL 

The model of color confinement is defined by the ex- 
pression for the free energy 

Wu{H) = c\k\V{[H^ + f{k)]"^ - f{k)}, (9) 

where c is a dimensionlcss constant, and f(k) > is 
an as yet undetermined function. For c = (2_D)^/^ this 
expression for Wk{H) satisfies the bound ([7]), as one sees 
from (iJ2 + /2)i/2 < \H\ + f. The boimd Q is saturated 
for small k if 



/(O) = 0. 
For example we shall later take 

/(fc)=m^-2|fc|, 



(10) 



(11) 



where to is a mass. 

One can easily verify that expression © also satisfies 
WkiO) = 0, as required for normalized probability p{A), 
and that the gluon propagator, defined below, is a posi- 
tive matrix, as required. 

By definition the "classical field" is the expectation- 
value 



ai{k,H)^{a';,{k))\H 



of the field 



^^(fc) = / d^x cos{k ■ x)Al{x) 



(12) 



(13) 



in the presence of the source H, and at momentum k. It 
is expressed in terms of Wk{H), by 



aiik^H)^ 



_ dWk{H) 



c\k\V 



h: 



[H^ + P{k)] 



and the gluon propagator by 



(15) 



P^(fc)J''-[g^ + /^(fc)]+g^gg 

Here P'^yik) = &jiv — k^k^/k"^ is the transverse projec- 
tor which appears because i?^ is identically transverse 
k^Hl = 0, and dHl/dH- = Pj^M^'^ 



The free energy Wk{H) vanishes at fc = for all i?. 



If Wk{H) were analytic in H in the limit fc ^ 0, this 
would imply that all derivatives of Wq{H) vanish and 
with it the gluon propagator at fc = 0, Djfi) = 0. 



However this disagrees with recent lattice data [1| which 
give a finite result, D{0) ^ 0, in Euclidean dimensions 



3 and 4. However the second (and higher) derivative 
W^'(iJ) is non-analytic in H in the limit fc — > when 
/(O) = 0, and the gluon propagator (fT5)) 



(17) 



does not necessarily vanish for fc — >■ 0. Indeed with 
/(fc) ~ TO-'^^^lfcl, we have 



D{0) = c/to 



D~2 



(18) 



which is finite. Thus our model, with /(fc) ~ TO^~^|fc| 
satisfies the bound d?]) and accords with lattice data 
for D = 3 and D = 4. In contradistinction to the 
previous treatment [§], the hypothesis that Wk{H) is 
analytic in H for fc — >■ is relaxed. Moreover, for 
/(fc) = m^~^+°'\k\^~", we get 



Dik) 



ck" 



(19) 



which fits the data in dimension £' = 2forQ; = l/5 p^ . 

If one fits the lattice data with this model, the data 
favor the value /(O) = in 13 = 2, 3, and 4 dimensions. 
With the value /(O) = 0, the model gives asymptotically 
at low fc 



Wk{H)^W^'{H)^c\k\\H\V. 



(20) 



This is linear in whereas normally one expects the 
free energy Wk{H) to be a power series in H, with leading 
term of order H^. 



III. MODEL QUANTUM EFFECTIVE ACTION 

As in statistical mechanics we define, for each momen- 
tum fc, the analog of the bulk magnetization 

dWk{H) 



M'{k,H) = 



(21) 



and make the Legendre transformation from Wk{H) to 

rkiMk)^MkH~Wk{H). (22) 

In fact the "magnetization" Mk{H), coincides in our 
gauge theory with the classical gluon field given in p4|) . 

ak{H) = Mk{H), (23) 

and the Legendre transformation (j22[) yields the quantum 
effective action of the gauge theory 



(24) 



More precisely this is the quantum effective action with 



(■^6) all variables set to except a^{k) for fixed fc. To find 



rfe(afc) we note that 
Tk{H) 



clklVH^ 



[i/2 + /2(fc)]l/2 



:\fik)v{ 



1 - 



-WkiH) 
/(fc) 



[H^ + P{k) 



TTj}' (25) 
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and from squaring (jl4p wc find 
/ 

[fl'2 + /2(fc)]l/2 - 



c\k\V 



1/2 



(26) 



which gives for the quantum effective action with all vari- 
ables except Cfc set to 

r,(a,) = c\k\f{k)v[l -[l- (^)'] "^}. (27) 



It is singular at 



\a{k)\ = c\k\V. 



(28) 
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